Stability of a vortex in a small trapped Bose-Einstein condensate 
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A second-order expansion of the Gross-Pitaevskii equation in the interaction parameter determines 
the thermodynamic critical angular velocity Qc for the creation of a vortex in a small axisymmetric 
condensate. Similarly, a second-order expansion of the Bogoliubov equations determines the (nega- 
tive) frequency uja of the anomalous mode. Although Sic ~ ~L^a through first order, the second-order 
contributions ensure that the absolute value \LOa\ is always smaller than the critical angular velocity 
Q.C- With increasing external rotation O, the dynamical instability of the condensate with a vor- 
tex disappears at SI* = \^a\-, whereas the vortex state becomes energetically stable at the larger 
value r^c. Both second-order contributions depend explicitly on the axial anisotropy of the trap. 
The appearance of a local minimum of the free energy for a vortex at the center determines the 
metastable angular velocity Q.m- A variational calculation yields Sim = \^a\ to first order (hence 
Q.m also coincides with the critical angular velocity Q.c to this order). Qualitatively, the scenario 
for the onset of stability in the weak-coupling limit is the same as that found in the strong-coupling 
(Thomas-Fermi) limit. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates 0-|^ in dilute alkali gases has stimulated great in- 
terest in the study of vortices in these systems. In the 
absence of experimental detection of a condensate con- 
taining a vortex, several authors have addressed the im- 
portant question of the stability of such a vortex [^|-|l6[ . 
The simplest criterion is energetic stability, determined 
by comparing the total energy of a condensate with and 
without a vortex, based on the time-independent Gross- 
Pitaevskii equation. Separately, numerical and analytical 
studies of the Bogoliubov equations for the excitations of 
a condensate with a vortex found a normal mode with 
negative energy and positive normalization, suggesting a 
dynamical instability |7|-|9|,p^. Moreover, in the limit of 
large particle number, there is a metastable regime, for 
which the free energy of the vortex condensate develops 
a local minimum at the center of the trap, stabilizing the 
vortex against small lateral displacements p 15 1. In the 
Thomas- Fermi limit, the onset of metastability coincides 
with |wa|, implying a close connection between the two 
phenomena |llf| . 

In the present work, we consider a small cylindrically 
symmetric condensate with a central vortex line along 
the z axis; this weak-coupling limit allows a perturba- 
tion expansion of the Gross-Pitaevskii (GP) equation. 
The additional energy of the condensate with a vortex 
can be compensated by rotating the whole condensate 
with an angular velocity Vic ^- The Bogoliubov equa- 
tions for the same system have an anomalous mode with 
negative frequency to a- To first order in the (small) in- 
teraction parameter, \ijja\ — ^c [p|p^- To understand 
the possible connection between these two types of insta- 
bilities, we here evaluate the second-order corrections to 



flc and uja- Furthermore, the metastable frequency Qm 
computed for the weak-coupling limit agrees to first order 
with the modulus of the anomalous mode, just as in the 
Thomas- Fermi limit. This result indicates that the in- 
stability with respect to microscopic oscillations and the 
onset metastability are closely related, for the system be- 
comes dynamically stable when rotated faster than the 
metastable frequency. 

The following section uses the GP equation for a con- 
densate with a central vortex to determine the critical ro- 
tation frequency ^c to second order. In the third section, 
a similar perturbation expansion of the Bogoliubov equa- 
tions for the anomalous mode yields the second-order 
correction to the corresponding eigenvalue of this mode 
LUa- The fourth section is devoted to determining the 
metastable frequency rirn variationally, and the conclu- 
sion discusses the implications of our results. 



II. CRITICAL ANGULAR VELOCITY 

Consider a condensate containing N particles in a har- 
monic trap with radial and axial frequencies lu± and 
LUz- The interparticle interactions are characterized by 
a positive s-wave scattering length a > 0. The starting 
point for the perturbation theory is the time- independent 
Gross-Pitaevskii equation [p7|jlq] 



with /i the chemical potential. Here, 
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is the Hamiltonian for the noninteracting condensate, 
expressed in dimensionless units (the radial and axial 



coordinates are scale d with t he radial an d axial oscil- 
lator lengths d± = y^h/Tnuj± and dz = y^h/muj^, and 
frequencies are scaled with the radial trap frequency 
hu;±). The condensate wave function ^ is normalized 
to unity, 7 = Na/dz is the small interaction parameter 
and A — uJz/uj± represents the geometrical asymmetry of 
the cylindrical trap. 

The normalized eigenfunctions for the noninteracting 
condensate are taken as a product of the one-dimensional 
axial oscillator state 1^9/(2) and the two-dimensional wave 
function for a g-fold quantized vortex x„+ q,„(r , (j)) at the 



center of the trap, so that [compare Eqs. (Al) and (A5) 
in the appendix] 



■ipn+q,n,l{r) = Xn+q,nir, (f>)(pi{z) 



(3) 



where ipqoo is the lowest energy state for a condensate 
with a given quantized circulation q. The noninteracting 
single-particle states have an energy e^^ ^ ^ ^ 2n + q + 

i + a + i)A. 

The thermodynamic criterion for stability of a rotat- 
ing condensate with a g-fold vortex is the vanishing of 
the free-energy difference A.Fq between the free energy 
of the condensate with and without the vortex. The free 
energy for a (7-fold vortex state in the frame rotating with 
angular velocity il is Fg = Eg — ^qN, where Eg is the 
energy of the condensate and q is the angular-momentum 
quantum number. Setting Eg — Fq = gives the critical 
rotation frequency 



Qr 



Eg — Eq 

Nq 



(4) 



To construct the perturbation theory, the chemical po- 
tential and the condensate wave function are expanded 
in the interaction parameter ^ « fi^*^' + 7^'^' -!-■■• and 

^ w vl/(o) + ^v]/(i) _| , The first-order correction to the 

chemical potential of the condensate is easily seen to be 
^(1) = 47r(*(")| |*(")p|*(o)), and the thermodynamic 
relation /i = dE/dN then gives the first-order correction 
to the condensate energy (note that 7 oc N) 



eW = lNfi<-'^ =. 27r7V(*W| |*(")|2 |*(0)) 



(5) 



For a q-fold vortex with ^'^'^^ — V'gOO: direct evaluation of 
the matrix element gives 



E, 



(1) 



N 






(2?)! 



(g!)2 229^/2^ 



(6) 



The first-order correction to the energy decreases with 
increasing q, and the corresponding critical rotation fre- 
quency 



nr 



1- 



E^^^ -E 
Nq 



(1) 



■7 



(7) 



increases with increasing q. Figure 1 shows the first-order 
approximation to the thermodynamic critical angular ve- 
locity as a function of the interaction parameter 7 for 



the first few values of q. Since the resulting Vie is always 
smallest for q = \, a singly quantized vortex will appear 
first, in which case fie ~ 1 — ^7/'\/27r. We consider only 
this case {q — 1) va. determining the second-order correc- 
tion to the critical angular velocity. 




FIG. 1. The critical rotation frequency in first-order 
perturbation theory for increasing quanta of circulation 
q = 1...6 (from left to right). The range of the interac- 
tion parameter 7 is extended to unphysical values for better 
illustration. 



The second-order correction to the energy of the con- 
densate can be found by writing '^^^'> as a sum over the 
unperturbed eigenfunctions of the appropriate symmetry 
(this form maintains the normalization to first order) 



*(i) 



/^ 



Cnl Xn+q,n ^21 



(8) 



where the prime on the sum indicates that it runs over 
all nonnegative values of {n, I}, omitting the single state 
{0, 0}. Substitution into the GP equation and use of the 
orthogonality of the unperturbed solutions determine the 
coefRcients 



27r 



Cnl = 



n + \l 
27r 



{lpn+q,n,2l\ IV'qOOp IV'qOo) 



n + Xl ' ««« 



(9) 



in terms of definite integrals discussed in the appendix. 
For example, the axial integration over the product of 
four harmonic-oscillator eigenfunctions yields 



dz Lp2l{z) <Pq{z) 



(-l)'./(20! 
^22' 



27r 



(10) 



with a similar but more complicated expression J"^ for 
the integral over the radial eigenfunctions (they involve 
Laguerre polynomials). 

The corresponding second-order correction to the con- 
densate energy follows immediately from the second- 
order chemical potential 



e: 



(2) 



N 



3^9 
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\{tpn+q,n,2l\ IV'qOOp IV'qOo)^ 



(0) 






= -87r 



^E 

71 1 



' (^i ^,"o?)' 
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(11) 



where /ig 



JO) 
tgOO 



9 



1 + ^A. The angular quantum 



number q is one for the singly quantized vortex or zero 
for the nonvortex ground state. 

In contrast to the first-order correction Eq in Eq. (|g), 
the second-order correction now depends explicitly on the 
asymmetry parameter A. The appendix describes a sys- 
tematic way of writing the matrix elements needed to 
carry out the infinite sums, and the difference between 
the second-order energy contributions for the vortex and 
nonvortex state Ai?^^^ yields the second-order correction 
to the critical rotation frequency 



^(2)/M^ 1 V' 1 l2 + 3n^-n^ {21) 
' ^ ' len^ n + Xl 22" 24'(/! 

nl 



(12) 



Consequently the second-order approximation to the 
thermodynamic critical rotation frequency (determined 
from the GP equation) is 



Qr 



1 



1 



2V27r 



7 + «(^)(A)7^ 



(13) 



l(2). 



Table I lists the second-order contribution 51c'' (A) for 
various asymmetry parameters A. This contribution is 
always positive, which counteracts the (negative) first- 
order contribution. For large A (a disk-shaped conden- 

sate), flc becomes constant (from the terms with I ~ 0), 
whereas for small A (a cigar-shaped condensate) , it grows 
like fie ex 1/A (from the terms with n — 0). The result- 
ing thermodynamic critical angular velocity for several A 
is shown in Fig. 2, illustrating the increasing importance 
of the second-order term for small A. For comparison, 
this figure also includes the critical angular velocity for 
a g = 1 vortex with only the first-order correction taken 
into account. 
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TABLE I. Second-order contributions for the critical ro- 



tation frequency Q 



(2) 



and the anomalous Bogoliubov mode 
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FIG. 2. The critical rotation frequency in second-order 
perturbation theory for decreasing asymmetry parameters 
A = 10, 2v^, 1, 0.1, 0.01 (from bottom to top). The solid line 
is Qc up to first order. Again, 7 is extended beyond the range 
of validity for illustration. 



III. ANOMALOUS MODE 

The preceding section dealt solely with the GP equa- 
tion for the condensate wave function ^ and the chemical 
potential fj,, comparing the energy for a condensate with 
and without a vortex. We now turn to the Bogoliubov 
equations [|l^ that describe the small-amplitude excita- 
tions of a condensate containing a singly quantized vor- 
tex: 

(Ho " M + 87r7 l^p) Uj - 4777 (^)^Wj = LUjUj , (14) 
-47r7 (**)^Uj + {Ho - f-i + 8nj |*p) Vj = -uj^Vj , (15) 

where Uj and Vj are the normal-mode amplitudes and Uj 
is the corresponding frequency. In zero order, the con- 
densate wave function 'i'^^' for a singly quantized vor- 
tex is given by -^loo ~ XioVo and the associated chem- 
ical potential is /i^"^ = 2 + ^A. For physical solutions, 
the normal-mode amplitudes have positive normalization 



JdV{\i 



\vj\ ) = 1. Here, we focus on the anomalous 



mode, which has a negative eigenvalue oja < 0. Hence the 
creation of quasiparticles in this mode lowers the energy 
relative to that of the condensate, implying a possible 
instability §-|,|l|,|§. 

A noninteracting condensate with a singly quantized 
vortex is unstable because the particles in the conden- 
sate with '0100 can make a transition to the nonrotating 
condensate with ipooo, giving up an energy 1 in our di- 
mensionless units. In the presence of interactions, how- 
ever, the description becomes more complicated, and the 
eigenfunctions of the Bogoliubov equation include both 
Ua and Va, even in zeroth order. At this level, the Bo- 
goliubov equations for the anomalous mode become 






(i/o-M^'^M «("' = -c^i"'«(°) 



(16) 

(17) 



where the preceding argument suggests that Ua oc 
"0000 = Xoofi^o characterizes the nonrotating vortex- free 
condensate, so that wi — — 1. Since the full Bogoliubov 
equation for Ua contains the coupling term (^)^ Va oc 
e^^'^fa, it is natural to assume that Va oc e~^"^, and the 
unperturbed state Va oc V020 = Xo2</5o has the correct 
energy to satisfy the remaining zero-order Bogoliubov 
equation. To ensure the proper normalization, we merely 
take 

u'^°'> = cosh6'xoo</'o and v^"^ = sinh6'xo2'Po , (18) 

where the parameter 6 can only be determined by includ- 
ing the higher-order terms in the perturbation expansion. 
The first-order terms in the Bogoliubov equations be- 
come 

-47r(vI/(«))\W = -(").« +c.W«i°), (19) 

-4.(vI/(«)*)M«)=-c.(°).W-c.W.r, (20) 

and we must also expand the parameter 9 « 6''-°-' +j6^^'' + 
■ ■ ■ . As a result, the zero-order functions in Eq. ( [l8| ) are 
evaluated with 9^^'. The first-order contributions then 
become 



M^i) = 0(1) sinh 0(°)xoo fo + Yl «»' ^"" ^2; , (21) 

nl 

!;(!) = 0(1) cosh0(°)xo2 <^o + y'bni Xn,„+2 ^21 ] (22) 



nl 



these expansions maintain the proper normalization to 
first order in 7. 

Multiply Eqs. (|l|) and (|^) by wi")* and wi"'*, re- 
spectively, and integrate. Straightforward manipulations 
yield the pair of equations 



-tanh0(o) I d?rxm{xiofxo2 



1 / tanh6i(o) 
2- 



a/2 



(23) 



A^W - C.W = 47r/o (2 1 d^^ |xo2nXioP 

-coth0(") I d^xhixlofxm 



1 /^3 COth^^O) 



n/2 



(24) 



where /q is given in Eq. (|lO|), and the radial integrals 
are evaluated with the results in the appendix. The sum 



of these equations determines IJ>^\ and comparison with 
Eq. d) yields 



cosh0(°) = %/2, and 
sinh6l(") = 1. 



(25) 



Correspondingly, the difference now determines the first- 
order correction to the anomalous frequency 



,(1) - 



2a/2^ 



(26) 



so that rif, = —i^a through first order. 

To determine the coefficients in the first-order expan- 
sions in Eqs. ^ and (||), project Eqs. (|l|) and d^) 
onto the appropriate unperturbed eigenfunctions. In this 
way, we find (see Appendix B) 



27r/, 



O-nl 



( / d^rxln{Xwfxo2 
2V2 /dVx^JxioPxoo) 



Ki = 



2Tr Ii 
n + XI 






(27) 



-2 



'^^''X;^,„+2(xto)^Xoo) 



27r/; / 



I \P) 7"2 _ O 7" 



n + \l 
The second-order Bogoliubov equations are 

(^(0)4^) +41)41) +4^).W) 

= (i7o-A^("))4^'-A^(^)..«-M(^)4"^ 
+87r [ |vI/(") l^^i) + (vl/(")*vl/(i) + M/W^/d)*) 40) 



(28) 



-47r[(*(0))2 4i)+2vl/(0)vl/(i)„(0)], 

-(4*')4^)+c.W4i)+4^).i°)) 

= (i/o-/.(°))4^)-4^««-M(^)«f 

-H87r [ |vI/(") p^i) + (vl/(")*vl/(i) + M/Wvi/d)*) i;(0)] 



(29) 



-47r[(*(°)*)2uii)+2*(°)**(i)*4°)]. 



(30) 



Here, Eqs. (||) and (g) give the condensate wave function 
*(i)^ and Eqs. (P), @, <^) and (||) give the correc- 
tions Ua and Ua ■ 

The remaining steps are essentially the same as in first 
order, yielding a pair of equations for ^(^) ±a;i ) . The pa- 
rameter 0(1) must be chosen so that the sum reproduces 
Eq. (11) for the chemical potential /i'^). The difference 
then gives the second-order correction to the anomalous 
frequency 



4^)(A)=87r^^' 



(Iif 



nl 



IX 



n2 \2 



[i3( jrA)' - siJifof MJiuf ~ 2(jri^o) 



_f jnO \2 Q j7il jnl A jnl jnl 



(31) 



and the total anomalous frequency through second order 
becomes 



1 1 



7 + 4'' (A) 7' 



(32) 



This is the main result of the present section, illustrated 
in Fig. 3. 
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FIG. 3. The absolute values of the anomalous mode fre- 
quency in second-order perturbation theory for decreasing 
asymmetry parameters A — 10, 2\/2, 1, 0.1, 0.01 (from bottom 
to top). The solid line is cja up to first order. 



the rotating frame is u!a{^) = i^a(O) + i^. With increas- 
ing external rotation fi, the condensate with a vortex 
becomes dynamically stable against microscopic oscilla- 
tions for n > U* = \u}a\- Our numerical results show 
that il* < ilc, so the onset of dynamical stability oc- 
curs before the singly quantized vortex in the condensate 
becomes energetically favorable. Since the sum flc + cua 
vanishes through order 7, it is necessary to include the 
second-order corrections to decide the relative value of 
fi* and fie- The difference is illustrated in Fig. 4 for two 
geometries, including the corresponding results for strong 
coupling (in the Thomas-Fermi limit) [0|, which shows 
that the sequence of the two types of stabilization is the 
same in both regimes. 
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FIG. 4. The difference A = f2c — \i^a\ (namely Qc + ^a) for 
the asymmetries A = 1 (diamonds) and A = 0.1 (crosses). In 
the weak-coupling limit, the difference is much smaller than 
in the Thomas- Fermi limit (note the logarithmic scale), but 
remains positive. 



For several different axial asymmetries. Table I com- 
pares the second-order correction V,c (A) to the criti- 
cal rotation frequency from the preceding section with 
the second-order correction wi ^(A) to the anomalous fre- 



,(2) 



quency. It can be seen that uja is always negative and 
smaller in absolute value than the corresponding correc- 
tion Qc ■ 

To understand the significance of this result, con- 
sider a normal mode of the q-fold quantized vortex with 
Uj{r, (j)) — exp[i(TOj-f g)(/)] Uj{r) and Vj{r, (j)) = exp[i(mj — 
q)(j)]vj{r) [|7|,pX|, where rrij is the angular momentum of 
the normal mode relative to that of the vortex in the 
condensate (namely, the corresponding perturbations in 
the density and velocity potential have the angular de- 
pendence ex expirrijcj)). In a frame rotating with angular 
velocity fi, the frequency of a given normal mode be- 
comes ujj{fl) = Wj(0) — rrijn. For the present case of a 
singly quantized vortex with q = 1, the anomalous mode 
has rria = — 1 , so that the shifted anomalous frequency in 



IV. METASTABLE FREQUENCY 

The third important frequency in this problem is the 
metastable angular velocity fim at which a slightly off- 
center vortex becomes trapped in a local minimum of 
the free energy at the center of the confining potential. 
Consider the corresponding energy functional 



E{^) = I dV [** {Hq - VtL^) * -h 27r7|*|*] 
We now use cartesian coordinates, with 



(33) 



H, = '- [-dl - dl + ix' +y') + \ {-dl + z')] (34) 

and Lz — —i{xdy — ydx) is the z-component of the 
angular-momentum operator. 

In order to evaluate the energy functional Eq. (^) , we 
construct a trial condensate wave function. This trial 
function is assumed to be unchanged along the axis of 



symmetry, which allows us to use the ground-state gaus- 
sian (po{z). In the weak-coupling limit, the radius of the 
vortex core is comparable with the radius of the conden- 
sate, so that the displacement Tq — (xq, j/o) of the vortex, 
the displacement ri — {xi, j/i) of the condensate, and the 
induced velocity of the condensate must all be included 
in the following trial function 



ipv{x,y,z) 



C 



[(x-a;o) + i(2/-yo)] 



X e 



-■^k^'-^iy' -^'yV^y^S^ ^2^^ J-{<^^^'+<^vV) 



(35) 



where C^^ = 1 + |ri — rop and the constants ctx and ay 
characterize the velocity components. 

In evaluating the integration in E{'>py), we retain all 
terms up to second order in the small parameters a^ and 
in the displacements. For convenience, we introduce new 
variables 5 = ri — yq = {6x,5y) and e = 2ri — ro = 
{ex,£y), which will turn out to be the normal modes of 
the system. The variational energy is 



Evi,r{ot,d,e' 



A 
2 



Q 



7 



2V27r 



+ — + ax {-5y + fifiy) + Uy {5 X - ^^x) 



The values 



n5 



— Uy 

= -Sx 



2n-'- + ^)s\ (36) 
2 V2^J 



Qe 



V 



(37) 
(38) 



minimize this expression with respect to the velocity pa- 
rameter OL = [axTay), and the energy then becomes di- 
agonal in the variables 5 and e, which thus represent the 
appropriate normal modes 






7 



/27r 



2n u^ 



i-n^ 



(39) 



The energy of the vortex-free ground state is E'q = 1 + 
\/2+j/{V2^). Thus the difference AF(0) = £;var(0)-£:o 
vanishes at the expected critical angular velocity Qc = 
1 — 7/(2\/27r), in agreement with the first-order result 
obtained from the Gross-Pitaevskii equation. Further- 
more, if S vanishes, then the vortex and condensate move 
rigidly, and the resulting dipole oscillation mode is stable 
throughout the range jfij < 1. For larger angular veloc- 
ity, the motion becomes unstable, as is familiar from the 
behavior of a classical particle in a parabolic potential. 

We next consider how the relative displacement S af- 
fects the variational energy Eq. (p9), which is clearly un- 
stable if O = 0. With increasing external rotation, how- 
ever, the relative displacement becomes metastable at a 
frequency 



i i^n 



7 



2V27r 



(40) 



when -Evar changes from a local maximum to a local min- 
imum for small relative displacements. To this (first) 
order in the interaction parameter 7, the value $!,„ coin- 
cides with the modulus of the anomalous mode frequency 
\uja\ and hence with the critical frequency flc- 

In order to illustrate the importance of the combined 
effect of the vortex and the condensate, we can study the 
behavior if only the vortex is displaced (namely ri — 0) 
or if the induced velocity is neglected {at = 0). In the 
former case, the energy (59^ is given by 



i^va.(ro)^2 + ^ + ^^ 



2 V2^ 



-n 



2n 



2 



(4 



vl) 



(41) 



and the metastable frequency is determined by requir- 
ing a positive coefhcient of the displacement contribu- 
tion [the second line in Eq. (E^)]. To first order in the 
interaction parameter, we find 



^*=l- 



7 



27r 



(42) 



which differs from the value fim found with the more 
general approach. In the second case, the metastable 
frequency is determined by the condition that the last 
line in Eq. (Bq) is positive, namely the determinant of 
the coefficients must be positive. This gives the same 
metastable frequency (^2|) as that found from Eq. (|4l|). 
Therefore, the metastable frequency has the same value 
as the absolute value of the anomalous mode (up to first 
order) only if the displacement of the condensate and the 
induced velocity are both taken into account. This indi- 
cates that the vortex becomes confined in a local central 
energy minimum at the same rotation frequency \uja\ for 
which the instability due to the anomalous mode disap- 
pears. The same scenario had been found in the strong- 
coupling limit pT| , and our results strengthen the idea of 
a common underlying phenomenon. 

Note that the metastable frequency found here 
[Eq. (pQ)] does not coincide with the metastable rotation 
frequency defined by Feder, Clark, and Schneider [[l5[ . 
Instead of requiring a local minimum in the free en- 
ergy, they identify the onset of metastability with the 
frequency for which the chemical potentials for a con- 
densate with and without a vortex are equal. In first- 
order perturbation theory, their criterion leads to the ex- 
pression (|4|) found by omitting the displacement of the 
condensate. 



V. CONCLUSIONS 

It has been shown that the critical rotation frequency 
U,c exceeds the modulus of the anomalous Bogoliubov 



mode uja when second-order corrections in the interac- 
tion parameter are included. This result agrees with nu- 
merical results for one particular trap geometry found 
by Feder, Clark, and Schneider [^. Furthermore, the 
sequence of stabilizing a singly quantized vortex through 
rotation is found to be the same as in the strong-coupling 
limit |1^: First, the vortex becomes stable against mi- 
croscopic oscillations at i7* = |wa|, and only at the higher 
rotation speed Oc does the vortex become energetically 
stable. The geometry dependence introduced through 
the second-order terms indicates that pancake geome- 
tries are more favorable for vortex detection in rotat- 
ing traps, since the stabilization frequencies are lower. 
This is consistent with the numerical results of Garcia- 
RipoU and Perez-Garcia [Q. For extreme cigar-shaped 
condensates, the critical frequency can even exceed the 
trap frequencies |15[| , a regime that is experimentally in- 
accessible. The angular velocity Q,m for the onset of 
metastability coincides with the modulus of the anoma- 
lous frequency tOa in first order, when the combined effect 
of the vortex and the condensate is taken into account. 
For VL > rim; an energy barrier stabilizes the vortex at 
the center; at the same frequency, the instability due 
to microscopic oscillations disappears. These features in 
the stability scenario also agree with the behavior in the 
Thomas- Fermi limit lllll . 
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APPENDIX A: NONINTERACTING 
EIGENSTATES 

The eigenstates for the noninteracting Bose conden- 
sate in a cylindrical trap can be classified with the quan- 
tum numbers of positive and negative circulation n_|_, n_ 
around the z axis and the axial harmonic-oscillator en- 
ergy quantum number I in the z direction [ pi[ . In par- 
ticular, the z-dependent parts of the eigenfunctions are 
simple harmonic oscillator eigenfunctions 



Vi{z) 



1 



■.Hiiz)e~-2' 



(Al) 



\/7ri/2 2'?! 

where the Hi{z) are the Hermite polynomials p2[ . 

In terms of the circular quanta, the normalized two- 
dimensional eigenfunctions are 



Xr 
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x^+y^ 



(A2) 



where al. = (aj, ± ia]^l\f2 = ^[x ± iy — {dx ± idy)] are 
the creation operators for right and left circular quanta, 
respectively. In terms of the new variables C. = x + iy, 
C* = a; — iy, these operators take the form 



t _C_ A t ^C_d_ 

"+ 2 dc* ' "- 2 ac 



(A3) 



The identity {^C " 9c)exp(-KC*) = exp(iCC*)(-5c) 
exp(— ^(^*) and its complex conjugate readily yield 
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Comparison with the standard formula for the associated 
Laguerre polynomials L^ |2^ yields 



Xn+q,n{^A) ^ (-1)' 



TT (n + g)! 



(A5) 



The complete normalized three-dimensional eigenfunc- 
tions of the noninteracting system are ipn+q.n,ii^) = 

Xn+qA'rA)'^l{z)- 



APPENDIX B: SECOND-ORDER MATRIX 
ELEMENTS 

The above expressions allow us to evaluate the neces- 
sary matrix elements, and the z-dependent part factor- 
izes out in the form 



dz Lp2i{z) ipo{zf 



71^/22' (2/)! 



dzH2iiz)e' 
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where only states with an even number of quanta con- 
tribute. The generating function for the Hermite poly- 
nomials 

-t^+2zt _ V^ Hn{z) 

readily shows that exp(— ii^) is the generating function 

for \/2'^^ {2iy. Ii , and a straightforward analysis gives the 
desired expression 
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For the radial part note that x„^„+g(r, (/)) — 
Xn-i-i3,n(^j (/)). Since the phase factor is the only complex 
part, it can be separated explicitly 



Xn,n+q{r,(t)) 

is 



X„+g,„(r), 



(B3) 



elements J"A. 



where x is real. The matrix 
/ (fr Xn+q,n Xifi Xj,o Xkfl involve a product of four of 
these eigenfunctions, where three refer to the conden- 
sate in the lowest energy state, namely n = (note that 
Lq = 1 for any a). Due to the angular integration, 



only p rodu cts with no net overall phase remain, 
of Eq. (|A5|) yields 
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where p = ^{q + i + j + k) is an integer because of the 
angular phase factors (note also that J"-^ is symmetric 
under interchange of its subscripts). As in the preced- 
ing example of the axial matrix elements, the generating 
function for the Laguerre polynomials |23] 
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exp 
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facilitates the radial integration. For example the func- 
tion 2(1 — t)/{2 — ty provides a generating function for 
(— l)"7r VrT+T J{\\, which is the matrix element with 
four singly quantized vortex eigenfunctions, and we find 
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A similar technique leads to all the other relevant radial 
integrals. 

The sums occurring in the second-order terms contain 
the square of the matrix elements (J"^ 7;)^ appropriately 
weighted with the energy denominator. For example, the 
second-order contribution to the chemical potential for 
the singly quantized vortex (obtained from the Gross- 
Pitaevskii equation) is 
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